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Abstract. We examine the behaviour of charged particles in homogeneous, 
constant and/or oscillating magnetic fields in the non-relativistic approximation. 
A special role of the geometric center of the particle trajectory is elucidated. In 
quantum case it becomes a 'fuzzy point' with non-commuting coordinates, an 
element of non-commutative geometry which enters into the traditional control 
problems. We show that its application extends beyond the usually considered 
time independent magnetic fields of the quantum Hall effect. Some simple cases 
of magnetic control by oscillating fields lead to the stability maps differing from 
the traditional Strutt diagram. 



PACS numbers: 03.65. Sq, 42.50.Dv 
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1. Introduction 

The present day quantum theories offer some visions of new, mathematically possi- 
ble reality which, until now, were not experimentally detected. Thus, the idea of 
supersymmetry looks for particles in superposed states of 'being fermion' and 'be- 
ing boson' and for the corresponding invariance group, with consequences for ele- 
mentary particles, strings, and quantum cosmology [3 [31 [31 IH [3 [3]. The other 
models propose to avoid the field singularities by assuming that the very struc- 
ture of the physical space forbids the exact point localizations. In some of them 
(Connes, Madore, Bellisard) the physical points have the non-commuting coordinates 
[71ll[Sl[inilIIl[Il[Il[Il[Il[Il[nilIl[IS]- According to Doplicher et al [H [H] , this 
might explain why the exact position measurements in the Heisenberg's microscope 
cannot produce the microscopic black holes, (even if the reality of the danger is an 
open problem). In the simplest model of the non-commutative plane (Madore [8]) the 
points {x, y) fulfill [x, y\ — in ^ 0, preventing as well the singularity creation. 

Though all these ideas are just free associations (in almost psychoanalytic sense), 
no better justified is the traditional hypothesis that our space is indeed a continuum of 
the 'exact points'. In fact, if adopting Wigner's observation about the "unreasonable 
power of mathematics" |22| . the mathematical models, if correct, must appear some- 
where in nature. So indeed happens, but with one amendment. The nature might 
'use' our models in its own way, without caring about the author's intentions. Thus, 
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while the problem of boson-fermion supersymmetry waits to be solved (still no trace 
of Higgs boson, no gravitinos, etc.), the same mathematical structure appeared in the 
supersymmetric quantum mechanics ( Witten [331 HI] ) helping to solve exactly a class 
of spectral problems (Duplij et al [B]). While there is still no sign of strings, branes nor 
extra dimensions, the analogous mathematical structures permit to understand better 
certain biological phenomena. In fact, the 'clones' of almost all unborn structures 
start invading the present day physics forming a little-great science in its own right 
(c.f. "big-bang in probet" etc. p5l[26]). 

Below, we shall discuss a similar status of the 'non-commutative points' which 
appear in the quantum Hall effect. We shall show that the idea works not only for 
the static fields but also for time dependent ones, even though the more fundamental 
hypothesis concerning the granular structure of the space itself is still waiting to be 
confirmed (or abandoned?). 



2. The charged particles in homogeneous magnetic fields: circular and 
drifting trajectories. 

We start from the known facts. The time independent magnetic fields B in an 
open domain of IR"^ can be described by a class of the vector potentials A(x) with 
B(x) = rotA = V X A(x). If B(x) is homogeneous, B = (0,0, i?) (for convenience, 
let B define the z-axis in IR'^) then one of natural choices of A(x) is: 



A(x) 



1, 



-B X X 



1 



-B 



X 





(1) 



interpretable as the vector potential created by a homogeneous current density on 
the cylindrical surface. The particle motion along the z-direction is then free; hence, 
we shall be interested only in the motion trajectories on the a;, y-plane. In the non- 
relativistic approximation, the Hamiltonian of a classical point particle of charge e 
and mass m is 



or equivalently: 



1 
2m 

1 
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[{P^ + + {Py - PXY 



(2) 



2^ ,r , {Plni)M. (3) 

where /3 = x = (x, y) and p = {px,py) are the generalized momenta, and 
— xpy — yPx- The well known shape of the motion trajectories is most easily 
derived from the Hamiltonian in form The first pair of the canonical eqs. defines 
the interrelation between the generalized and kinetic momenta mv 
dx dHjj 1 
dt „. 

-{Py — PX) Py — mVy + f3x 



dy 
dt 



dpx 



m 
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and the second pair yields the proper dynamical equations: 

dpx _ dHp _ P_ 
dt dx m 

dpy _ dHji /3 

dt dy 



-{Py- Px) 

-iPx + Py) 



(5) 
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Both (H])-® imply immediately the existence of two conservative quantities: 
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(6) 
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(8) 



so, each charged particle just rotates around a fixed center X = (X, Y) with a constant 
(cyclotron) frequency a; = ^ = — : 



x{t) - X 




yit) - Y 





(9) 



cos(a;i) sin(a;i) a::(0) — X 
— sin(a;i) cos(wt) y{{)) — Y 

The conservative quantities X, Y represent the hidden symmetries of the system 
[TUl [T51 [13] • The expression for X, Y in terms of the generalized momenta might 
look peculiar but it turns natural in terms of the velocities (kinetic momenta) (|4]). 
It shows that the radius p of each rotating trajectory depends just on its (constant) 

velocity scalar w = |v| = y^(i)^~+~i)2): 

-,2 



p' = {x-XY + {y-YY = 



(10) 



Yet, to find the system response to the external forces, the most convenient are the 
expressions dSHll) in terms oipx^Py Indeed, if the circulating charge (IS])-® is affected 
by an additional potential V^(x), the Hamiltonian becomes H = Hp + V"(x) and since 
the motion center is conserved by i?^, the canonical eqs. for X, Y are reduced to: 



/3' ^ '^V 2/3 9y 2/3^ 



(11) 



(12) 



dt ^ ' ^ 2^ dx 2p 
In the simplest case, if the trajectory ([9]) is affected by a constant force F, the well 




Figure 1: The three generically similar phenomena, a) The rectilinear drift along 
the axis y under a constant force F in the x direction, b) The circular motion 
trapped by a repulsive oscillator, c) The precession of a gyroscope. 
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known though counter-intuitive effect is that the rotation center starts to drift in the 
direction orthogonal to F. Curiously, the exact solutions of exist also for the 
elastic potentials V^(x) — ax^. For a < 0, \a\ < they illustrate the surprising 
phenomenon of a charged particle trapped by the repulsive center (Fig. [1]). 

A less elementary form of the same effect, in presence of crossed electric and 
magnetic fields is discussed in [27i, where the KAM theorem [28| is applied to show 
that the existence of an additional repellent obstacle (in form of a disk) can interrupt 
the rectilinear drifting, trapping the charged particle which instead of being rejected, 
'obsessively returns' to the repelling obstacle. Below, we shall not pretend to deepen 
this line of thought, but instead, we shall focus attention on the corresponding 
quantum systems. 



3. Quantum case: the "fuzzy points". 

"Pure quantum states are objective but not reaV 

Hans Primas, Chemistry, Quantum Mechanics and 
Reductionism., Springer Vcrlag (1983) p. 103. 



The quantum equivalents of p ll2p are so widely studied, that we fix attention 
just on one particular aspect. We shall use the same symbol Hp to denote the 
quantum Hamiltonian ©-(O, though now Xk^Pi will mean the quantum observables, 
[x/.,Pj] = iSki- By commuting the Hamiltonian ^ with Xk,Pi one obtains the same 
two pairs of equations ^ and ([5]), though now they concern the Heisenberg's operators 
x{t) , Px{t) , y{t) , py{t) . In a complete analogy with the classical case, the formulae (|4])- 
([5]) show the existence of a pair of conservative observables X, Y given by . 

The works on quantum systems ©-© dedicate a lot of care to the non- 
commuting kinetic momenta [7J 129) . It neither escapes attention that the same 
phenomenon affects the coordinates of some abstract 'space localizations' such as the 
rotation center of the Heisenberg's trajectory (c.f. Avron et al. [3D], Dodonov et al. 

my- 

y-\v-Yp^^ (14) 

While the instantaneous particle coordinates x{t),y(t) at any fixed time moment t 
are still the commuting observables, this is no longer true for the rotation center 
X — {X,Y) which becomes a "fuzzy point" ]8^: 

[X,Y]^~^^-— (15) 

Note that, the "fuzzy localization" is an abstract concept (resembling a "grin without 
the cat" of Lewis Carroll [32] )• There is simply nothing there, commuting or not. Yet, 
by paraphrasing Primas |33| . the fuzzy center (jl3p - (|14p is "nof real but is objective^\ 
Another curious aspect, resembling some open gravity problems |341 1351 136] 
concerns the 'surface' of the 2-dimensional orbits (|9I10I) . Indeed: 



TT 



irp =^(x-X) = — 



(Px + Py) + (Py - 



;Hfi (16) 
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Hence, the surface of the orbit is not only conserved but also quantized: it is 
proportional to the Hamiltonian and cannot change continuously (an equivalent 
observation see, e.g. j57]). 

As it seems, this is not the first time when the magnetic fields provide an imitation 
of still unchecked theories. The suggestive ideas of supersymmetry are still not verified 
in particle physics. Yet, due to the anomalous relation between the spin and orbital 
magnetic moments of the electron, its energy levels in a static, homogeneous magnetic 
field reproduce the supersymmetric spectrum. So, many authors conclude that: "the 
supersymmetry exists in nature" [351 IMl HSl SI] ■ Quite similarly, there is no evidence 
that the physical particles move in non-commuting spaces, but the properties of 
the rotation center (jl5|) might encourage the statement that "the non-commutative 
positions exist in nature" . 

It might be thus interesting to notice that the phenomenon is not limited to the 
circular trajectories in the static magnetic fields. It appears as well in more general 
physical scenarios, including the time dependent oscillator potentials and magnetic 



4. Time dependent oscillators: the non-circular loops 

The evolution problem for variable oscillator potentials has already a notable past. 
The first systematic studies were presented in 1967-69 by Lewis and Riesenfeld [42ll43j . 
then by Malkin and Man'ko [IJ, inspiring an ample research on the state evolution. 
The group theoretical approach to the Baker-Campbell-Hausdorff (BCH) problem 
PSI ITTl H5] for quadratic Hamiltonians was developed in 1974 by Gilmore 05] (a 
generalization see Zhang et al [SD]). In 1976 Yuen notices the possibility of solving the 
evolution equation for the squeezed photon states [5T]. For massive charged particles 
the operations induced by variable fields were studied by Ma and Rhodes [52], Royer 
[53] . Brown and Carson [54], Combescure [55], Wolf [56] and other authors; for an 
ample review c.f. Dodonov [57]. The links with the coherent states are notable 
[58l[59l[60ll61]. In the parallel research there was hardly any constant in the generalized 
oscillator or magnetic Hamiltonians which would not be replaced by its time dependent 
analogue, starting from the time dependent mass {c.f. the Hamiltonian of Caldirola- 
Kanai [BHIMIISI]) up to the variable dielectric or magnetic permeability [551 [551 1571 IB5] . 

The possibility of the closed non-circular trajectories for the time dependent 
oscillator Hamiltonians was noticed already in Malkin and Man'ko [H] . The extremely 
simple cases of the closed evolution caused by the sudden (5(i)-shocks of the oscillator 
potentials were described in [69l [TOl [TlJ [72]. Thus, for the Schrodinger's particle 
in 1 space dimension the evolution loops can be produced by sequences of oscillator 
pulses and free evolution intervals corresponding to the elementary cases of the BCH 
formula, e.g. 



fields. 
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3± ^ 
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(17) 



6 terms 



(with h = 1), illustrated by the evolution diagramme: 
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where the vertices symboUze the shocks of the elastic potential (with the corresponding 
numbers, meaning the pulse amplitudes) and the sides correspond to the 'rest 
intervals' of the free evolution (for simplicity we accept the particle mass m = 1). 
The equivalence sign = in (jl7p means the operator proportionality, i.e., U = 
1 <J4> [/ = e^"^!, ip & R. The many vertex analogues of (|17l) as well as the 
more general 'kicked systems' and their non-singular equivalents were described in 

Some incomplete versions of the loop process may prove of interest. In fact, 
whenever an evolution loop contains a i5-pulse of the attractive oscillator potential, as 
in (fTTl) , the whole rest of the process must imitate the effect of a repulsive elastic pulse 
(which cannot be straightforwardly achieved by the magnetic fields). Furthermore, if 

_■ eI 

any evolution loop contains an interval of free evolution e 2 (r > 0), it means that 
the whole rest must be equivalent to its inverse, thus suggesting the techniques of 
reverting the free propagation. The importance of the closed dynamical processes for 
more general control operations was recognized in [69l [72l [74j [75l [78] . 

All these models depend on some idealizations. Thus, the repulsive oscillator 
kicks cannot be straightforwardly simulated by the magnetic fields (in fact, even 
the attractive ones can hardly be engineered!). Moreover, in order to describe 
the time dependent potentials in some little but macroscopic areas as, e.g., an ion 
trap, one typically uses the non-relativistic laboratory approximation, disregarding 
the little delays needed to propagate the potential inputs all over the trap surfaces 
(telegraphist's equations) or in its interior. Since these delays in typical experiments 
are insignificant, the lab approximation works very well and indeed, is implicit in all 
papers postulating the time dependent external parameters {e.g. [42l EH [TlJ [TOl l82j 

iHES]). 

Once this approximation is adopted, one can see that the fuzzy centers are not 
restricted to the orthodox Hall effect with a fixed magnetic background. They arise, 
with comparable consequences, for arbitrary evolution loops generated by the static 
or time dependent quadratic Hamiltonians. 

5. Loops of quadratic Hamiltonians: stability and drifting. 

To show this, we shall consider the hermitian Hamiltonians H(t) defined by time 

2s 

dependent quadratic forms H[t) = ^ hij{t)qiqj (with hij(t) = hji{t) G R) of a 

complete set gi, • • • , q2s of any number of the canonical observables, Xi,pi {i — 1, s) 
in a certain Hilbert space Jf; [xi,Xj] = \pi,Pj] = 0, = iSij. An agreeable 

property of the quadratic Hamiltonians is that even though H{t) are unbounded, but 
if the coefficients hij{t) are non-singular and piece- wise continuous, the corresponding 
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unitary evolution operators U{t, t) are well defined by the operator equations: 

j^U{t, t) = -I H{t)U{t, t) ^ ^U{t, t) = I U{t, t)H{t), (18) 

with U{t,t) = 1, and: 

U{t,T)U{T,a) =U{t,a), (19) 

(deeper results and amendments, c.f. B. Simon |83j, Hagedorn et al ^Mj)- As well 
known, for quadratic H{t), the time dependent Heisenberg's observables Xi{t),pi{t) 
are linear combinations of the initial Xi,pi. For convenience, we shall denote by q the 
vector-column of 2s (dimensionless) observables: 





qi 




X 


q = 














p 



(20) 



One therefore has: 

q{t) = U{t,OyqU{t,0)=uit)q (21) 

where u{t) — u{t, 0) is a dimensionless 2s x 2s evolution matrix yielding simultaneously 
the classical and quantum trajectories. 

If q generates all observables in J^, then not only the evolution operator U{t) 
determines the matrix u{t), but inversely, the matrix (|2ip determines the unitary U{t) 
up to a numerical phase factor. Indeed, should U and U' be two unitary operators 
generating the same transformation of the canonical var iables q, i.e. WqU = U'^'qU', 
then U'U^q = qll'W, and so, U'U^ would commute with all canonical variables x, p 
and their functions. However, if the algebra spanned by the canonical observables q 
is irreducible in J^, then any operator which commutes with all of them must be just 
a c-number. Since U and U' are unitary, this number can only be a phase factor, 
i.e.: U'W = e"^ implying U' = e^^C/. The concrete value of ip is relevant for the 
transformation of the state vectors but not of the Heisenberg's variables, so we shall 
simply denote: 

U' = e'^PU => U' = U (22) 

This equivalence turns quite essential in quantum control, as it basically permits to 
deduce the form of U from the transformation of the canonical observables. Moreover, 
since for the quadratic H{t) the evolution matrix in quantum and classical cases is 
the same, the properties of the evolution operators U{t) can be read as well from the 
classical motion trajectories. This allows numerical but non-perturbative solutions of 
the continuous BCH problem [Ml |49l |86l |87] for the evolution (fT8|) . In particular, an 
evolution loop occurs if after a certain time interval (for convenience, let it be [0, T]) 
all canonical (classical and/or quantum) variables q return to their initial values: 

q(r) = c/(r,o)tqc/(r,o) = q ^ ^^(t) = i, (23) 

so that, in the sense of (P^ : 

[/(T,0) = e^'^l = 1, (v?GR), (24) 

implying also a loop of all other observables A{t) which do not depend explicitly on 
time in the Schrodinger's frame, i.e., A{T) = [/(T, 0)U[/(r, 0) = A. The value of 
ip £ R {i.e., the geometric phase for the loop process) though interesting in itself 
[88l[72], does not affect the results of our present argument. 
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Below, we shall consider a Hamiltonian H{t) varying periodically, i.e., H{t + T) — 
H{t). We denote for simplicity U{t,0) = U{T + t,T) — U{t). We assume moreover 
that the evolution produces a loop in [0, T] and in the subsequent periodicity intervals. 
Then, any Heisenberg's observable A{t) = Uit^ AU{t) in any periodicity interval 
[r, T + T] admits the time average 

T + T 

A=^ j A{t)dt, (25) 

r 

independent on t, defining a certain global characteristic of the loop process. (Indeed, 
the change of r in ()25|) means just that the same contributions A{t)dt are summed up 
along the same closed cycle, rearranging only the summation order jj. 

In particular, applying ([25)1 for the time dependent variables Xj{t), one obtains 
the coordinates Xj of the 'loop center'. For H{t) quadratic they are as well linear in 
the initial qi, ...,q2s- 

X = X ^ Xj ^Xj ^ Ujiqi H h 2s '?2s = 1, • • • , s), (26) 

generalizing the already described rotation centers (IT3)) - (IT4l) in 2D. A more general 
concept might be also useful. 

Definition. For the quadratic, periodic Hamiltonians, even if the evolution in 
[0, T] does not close to a loop, one might apply I^E^ defining the Floquet point 
X(r) — y JJ^'^ x(t)dt (though, in general, X(r) will depend on r). 

We shall show now, that the behavior of the loop affected by an additional force 
described in Sees. [5][3]is a typical phenomenon for all loop processes generated by any 
quadratic Hamiltonians. 

Proposition 1 Suppose, a time-periodic Hamiltonian H{t) = H{t + T), quadra- 
tic in the canonical ohservahles q, generates a loop \24^ in its periodicity intervals 
[nT,{n + l)r], n = 0,1,.... Then the precession or stability of the loop under an 
additional, constant force F depends on the loop center X = (Xi,---,X„). // its 
coordinates vanish, then the loop can change its shape hut remains closed and stable. 
If the coordinates commute, then the center is stable even if the trajectory is not. 
However, ifK. is a "fuzzy point" with non-commuting coordinates, [Xk,Xi] — iKki ^ 0, 
then the loop will show a drift with a constant velocity ^ Vk = HkiFi (summation 
convention) in the direction orthogonal to F . 

Proof. Suppose the perturbing force is F — {Fi,...,Fs). The perturbed 
Hamiltonians then read: 

H{t) = H{t) - Fx, (27) 

yielding the modified evolution operator U {t, 0) in the customary form: 

U{t,Qi) ^U{t,Q)W{t) (28) 

in which U{t,Q) is the evolution operator of the unperturbed loop ([M)) . while W{t) 
obeys the evolution equation in the interaction frame: 

^ = iF^{t)W{t); W{Q) = 1 (29) 

I This does not mean that A is a constant of motion in a conventional sense. As an integral over 
the entire time interval [t, T + t], the observable A is not 'local in time'. Yet, if the little evolution 
steps of A{t) in I I25II obey their (different) instantaneous tiamiltonians, then II25I I stays unchanged. 
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where 

x(t) = C/(t,0)^x[/(t,0) (30) 

are the time-dependent Heisenberg's observables defined by the unspoiled loop 
evolution (fT8|) . (f24)) . Since the commutators [xkit), xi{t')] for any t,t' G R are numbers, 
the operators W{t) are given by the simplest case of the BCH formulae [85], [86l |87] 
in which only the integral of x(t) matters: 

W{t) ^ eM'^e'^J'o'''-''^'^'' , (31) 

and x{t) is a real, c-number phase for each t G R. So, 

U{T, 0) t/(r, 0)W{T) = e'^^e*^^^, x e R, (32) 

implying: 

C7(T, 0)^ C7(r, 0) = e~*^^^FXe*^^^ FX, (33) 
but simultaneously 

UiT, 0)^YU{T, 0) = e-*^^^Fe''^^^ = F - [FX, Y]=Y + oT, (34) 

for any Y = niXi, where D = niXijFj (summation convention). Hence, the (constant) 
external force F cannot change the center coordinate FX, but if [FX, y] ^ it can 
shift the combination Y of the remaining ones. If the loop is n-times affected by the 
same F, then FX is still unchanged, but Y performs a cumulative drift: 

U{nT, 0)^YU{nT, 0) = W {T)^'^YW (jy = e-^^^^Fe™^-^^ = Y + nTv (35) 

with the same constant velocity V.m 

Note. Above, it is not assumed that T is the smallest period of H{t), but only 
that it is a common period of H{t) and of the loop phenomenon. In fact, a typical 
situation is that the loop occurs after several periods of H{t) (compare Sec. [71). 

Our results permit to classify the loops according to the properties of X: 

• A loop generated by a quadratic Hamiltonian for s > 2 can have a 'fuzzy center' 
X with non-trivial [XfejA";] = inki ^ 0. Then under the influence of a constant 
external force F the loop can be broken; the former motion center will drift with 
a constant velocity in a direction orthogonal to the applied force. 

•• If [Xfe,X;] = 0, then the loop center X, no matter the external force, must return 
to its initial value. Various subcases are still possible. 

i If X does not vanish identically, the loop still can be broken. The motion 
center won't drift but the trajectory can precess around the fixed center 
X. The detailed classification of the loops in this subcase is still an open 
problem. 

ii If X identically vanishes, then W{T) in (|32]) is just a phase factor and the 
loop is stable: it can be deformed but it won't be broken, neither will it 
precess under the influence of F. 

The sense of this classification turns obvious if one compares the traditional harmonic 
oscillator with the "magnetic oscillator" Hp (/3 =constant). Both admit circular 
orbits, however, the general elliptic orbit of the oscillator affected by an external force 
F gets simply displaced in the direction of F (where it remains stationary), while the 
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orbit of Hp starts drifting in the direction orthogonal to F. The key to this difference 
is the distinct nature of the classical/quantum center {X,Y) of both motions. While 
for the 2D oscillator the center X is exactly zero, the same center for Hp is a fuzzy 
point P^ - (|14p implying the drifting trajectory. 

In physical terms, the difference between both cases has some historical key. In 
fact, by reading (|3]) inversely: Hose = Hp + ^M^, one obtains a modern equivalent 
of an old idea: the description of the elliptic orbit of i?osc as the superposition of two 
circular motions, for Hp and M^, i.e. the Ptolemean picture of the oscillator trajectory. 
So, in a sense, "the epicycles are more stable than cycles" . 

6. Landau fields: manipulating tlie fuzzy center 

As already known, the magnetic operations can produce special effects such as the 
rigid displacement, squeezing, and the distorted free evolution (TTJ |^ I^H [31] . Below, 
we shall be specially interested in the homogeneous magnetic fields in a fixed direction 
(e.g. of z-axis), imitating the variable oscillator potentials in 2D [77l [3T1 [70l [7T1 [89] . 
The simplest physical conditions to approximate such fields arise in space domains 
surrounded by time dependent currents, e.g., by solenoids of various forms (Fig. [2^,b). 




Figure 2: a) Cylindrical geometry, b) Landau geometry. 



If the surface currents do not depend on z, the motion in the x, y-plane decouples 
leading to some typical vector potentials in 2D, e.g.: 

B 



cylindric : 



A = 



-y 



Landau 



A = B 



-y 





(36) 



If B is static, both expressions (|36l) are exact and gauge equivalent. Yet, for variable 
B = Bit) they still offer a good laboratory approximation, generally used in the ion 
trap description, though neglecting the retarded field effects in not too huge operation 
areatH. Below, we shall check the loop behaviour for the time dependent quadratic 
Hamiltonian of the Landau's case: 

2 



^ ^ 2m 



eB{t) . 2 
Px + ^^y] +pI 



(37) 



§ Differently than for the static B, the cyUndrical and Landau potentials in I I36II for variable B{t) 
are no longer gauge equivalent, since they produce different electric fields, corresponding to distinct 
geometries of the distant sources; c.f. Dodonov et al. | 31| . 
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In the diniensionless variables with t — >■ i t and /? = 
unit), it reads: 



e_B(t)T 
2mc 



Hit) = i 

The canonical equations 
dx , 



{p^+2mv) +pI 



(where t represent a time 



(38) 



2/3y) 



dpx 
dt 







dy_ 

dt 
dpy 
dt 



Py 



^ = -2/3(p. + 2/3y). 



(39 a) 
(396) 



despite their apparent simplicity, require a computer study. We thus opted to 
approximate an arbitrary /3{t) by a step function of 2n constant fields /32, ■ ■ ■ : P2n 
in the subsequent 2n time intervals ti], i — 1, ■ ■ ■ ,2n {0 ^ to < ■ ■ ■ < t2n — T). 
While such steps are never exact (to switch on or switch off an electromagnetic field 
takes at least 100 attosecond, [90], [9TJ[92]), our aim here is just to check the principal 
part of the behaviour. To simplify still further, we fixed the Pi and Ati = ti — ti^i 
to convert each step into a '7r-pulse' generating a certain semi-circle in the x, y plane, 
with PiAti = ±^ (« = 1, 2, • • •). The semicircular fragments of the canonical trajectory 
correspond to the evolution matrices 



Ui = u{ti,ti^i) 





_ j_ 
1 



Hi 




(40) 



leading to the final result of 2n operations: 



u(T) = U2n • • • wi = 1 + r 



2n-l 

E 

i=l 



1 



1 



(41) 



b) 




Figure 3: a) In black: the classical trajectory illustrating the loop (|4ip with 
ft = f , Ati = 3; ft = f , At2 = 2; ft = tt, At3 = i; ft = f , At4 = |. 
The "fuzzy center" is (X, F) = [x + j^Py, --^Px + j^Py) ■ In blue: the force 
F = (0, 1) causes the drift of the loop center, b) The loop for Pi — j, Ati = 2-k; 
ft = 1, At2 = f ; ft = -J, At3 = 2tv- ft = -1, At4, = f . The motion center 
(X, y) = (a; + Ipy, -Py) is operational but no fuzzy. Green: affected by the force 
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with the loop condition u{T) = 1 4^ T = 0. 

The simple form of the semicircular operations makes possible not only to 
draw the trajectories but also to calculate exactly their centers X — {X, Y) and 
commutators [X, Y] . As an elementary example, we did it for the 4 step loops 
obtaining: 



+ 



1 



+ 



1 



(42) 



L/3ii/3ii /32I/32I km MM 

Within the loop condition F = 0, the commutator (j42|) admits still various control 
options, producing the "fuzzy centers" with distinct drifting capacities (see Fig[3K)- 
However, [X, Y] may also vanish for some values of /3i. This happens e.g. for (3^ = — /3i, 
Pi = —(32 when the quantum center is operational, but not fuzzy (see Fig. |3l3), the 
case which cannot occur for static fields. 

From the relativistic point of view the Hamiltonians p7p are inexact, since the 
retarded effects are missing. Notice, however, that the non-relativistic Landau's 
potential p6l) with B = B{t) is the limit for -i — > of an exact, relativistic expression 
in form of a finite difference: 

A{t,^) . ... ..^ lG(i~f)-G(i + f) 



A(t,x) = 







A(i,x) = -^ 



-G'{t)y 



(43) 



of two plane Landau pulses propagating in opposite directions, where G{t) is a 
continuous, real function with a bounded, piece- wise continuous derivative G'(t) — 
B{t), modeling the time dependent magnetic fields for i — >■ 0. So, the non-relativistic 
Hamiltonian (l37t with the time-dependent B{t) is indeed the 1"* step of the Einstein- 
Infeld- Hoffman (EIH) approximation commonly used in General Relativity to describe 
the slow motions in limited space domains |93j . In the laboratory of the size of « Im 
the y-dependent delays don't exceed St = = | • 10~*sec and the approximation is 
almost perfect (errors invisible in our Fig. [3^,b). We thus conclude that in not too 
huge magnetic traps the "fuzzy" centers of the Hamiltonian (|37p describe correctly 
the principal part of the drifting phenomenon. 



7. Cylindrical geometry: the resistant loops 



In turn, we shall examine some curious phenomena induced by softly pulsating fields 
of cylindrical geometry. The non-relativistic Hamiltonian becomes: 



H{t) = ^ 



eB{t) 
2c 



I 2rnc ) ^ 



(44) 



Differently than the 'homogeneous Schrodinger's equation' (c.f. [Il]), the eq. (|44| 
represents a 2D 'Aristotelian world', whose symmetry center x = y = is distinguished 
by the circular electric fields (c.f. also [SI])- In principle, B(t) might be arbitrary, but 
below, we shall be most interested in the harmonic and biharmonic fields: 



B{t) = Bo + Bi sin(wt). 



(45) 



B{i) = Bi sin(wt) -I- B2 sin(2a;t), (46) 

of period T — Following the commonly applied approximation [331 [Ml [57] we 
adopt the semiclassical picture (Thomson rather than Compton [95]). Our design is 
quite unsophisticated comparing with the time dependent mass [621 163j or the variable 
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material constants [681 167) . However, what precisely happens with the microparticles 
in this simple scenario? 

To describe the ample classes of similar evolution processes it is practical to 
introduce the dimensionless time, t ^t — 4^, field B = -^B, and canonical 
variables [a;,^^;] = = The field oscillations now become: 

(3{t) = ^0 + /3i sin(27rt) + ^2 sin(47ri), (47) 

where ,92 = corresponds to the harmonic and /3o = to the biharmonic cases (H5|) 
and p6)) respectively. The rescaled Hamiltonian is 



H[t) = \ 



+ /3(t)2x2) - /3(t)M, , 



(48) 



Hose Hj^ot 

with -ffosc representing the time dependent "magnetic oscillator" and Mz the rotation 
generator. Since Hosc{t) and Mz commute, U{t) factorizes into two commuting unitary 
operators, U{t) = Uosc{t) Urot{t), where: 

dUoscit) 



dt 



— I- Hose 

n 



(t)c/osc(i); 



C/osc(0) = 1, (49) 



while 



(50) 

produces just the rotations r(t) between the canonical pairs x,Px and y,Py The 
canonical transformation defining the evolution matrix u(t) can be split into two steps: 
q — >■ qosc(i) — >■ q(i), implemented by Uosc{t) and J7rot(i) respectively. The operation 
Uoscit) is reducible, affecting separately both canonical pairs x,px and y,Py which 
evolve simultaneously according to the same 2x2 matrix further denoted b{t), i.e.: 



X 
Px 



Uoscit) ^b{t) 



X 
Px 



(51) 



idem for y,Py By differentiating both sides of (1511) in agreement with (j49p and using 
the canonical commutation rules, one sees that b{t) is determined by the differential 
matrix equation: 



db 
dt 



A{t)b{t), A{t) 





-Pit? 



6(0) = 1 



(52) 



which is at the bottom of all quantum control problems for time dependent oscillator 
Hamiltonians [42l |43l HH [53l [TOl Ell [96] . Once having 6(t), one immediately constructs 
Uosc(i), as the simple pair of two 6(t)-cells. In turn, multiplying Uosc(i) by the 4x4 
matrix r(i) rotating by 7(i) — /3(i')c?t' around the z-axis, one obtains the complete 
4x4 evolution matrix 



uit) = rit) 



bit) 



bit) 



(53) 



Here the role of bit) is specially relevant for the periodically repeated field 
patterns, since the one-period-evolution step 6(T) decides about the bounded (stable) 
character of the motion, or its capacity to produce the parametric resonance 
[9711981 El]. 

In our case (|47p . T — 1. As representing the canonical evolution, 6(i) are 
simplectic, including the Floquet matrix 6(1), which permits to classify the motions. 
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Since Det[6(l)] = 1, its eigenvalues depend on just single (real) trace invariant. The 
characteristic equation 



D{X) = Det(A- &(1)) = - ATr6(l) + 1 = (54) 



has two non-vanishing roots 




with A+A_ = 1, distinguishing three possible types of motion: 

I. Stability area. If |Tr6(l)| < 2, then A+, A_ are two different complex eigenvalues 
with |A+| = |A_| = 1 (phase factors) of the form A+ = e*"^, A_ = e~"fi,ipe R. The 
6(1) shows an oscillating behaviour of b{n) = 6(1)", n e N defining the bounded 
trajectories. 

II. The threshold (separatrix) is characterized by |Tr6(l)| = 2. The Floquet matrix 
6(1) here has two coinciding eigenvalues A+ = A_ = ±1. If 6(1) is diagonaliz- 
able then once again b(t) trajectories are bounded, but if not, then they can show 
a weak parametric resonance growing in arithmetic but not geometric progression. 

III. If |Tr6(l)| > 2, then 6(1) has a pair of real eigenvalues, A+, A_ 7^ 0, A_ = j-^, of 
which at least one has the absolute value > 1. The trajectories show a strongly 
resonant behavior of the squeezing type. 

In the best known case of ion traps with the sinusoidally oscillating elastic force 
the solutions b{t) are expressed in terms of the Mathieu functions and the resonance 
borders form the well known Strutt diagramme (see e.g. |100| I96j). Due to its familiar 
shape, most of studies illustrating the effects of the time dependent elastic potentials 
stick to the Mathieu scheme (see e.g. Paul's trap the quantum tomography 

[55] . etc.). The need of some wider stability designs was pointed out by Glauber 
(c.f. the statement in Baseia et al |101| ). Indeed, some distinct stability cases were 
considered for finite dimensional state spaces [1021 11031 1104) . or else, for the rotating 
magnetic |105| 11061 1107] or electric fields |108[ 1109] . However, this does not include 
the magnetic oscillators Hose (i48]) modulated by f3{t)'^. The difference seems modest 
but the consequences are not. 

7.1. Harmonic loops: an anomalous resistance? 

A systematic study of the harmonic magnetic case was undertaken by F. Delgado 
[5^1110) by scanning the stability thresholds for the time dependent magnetic oscillator 
(|48|l with the harmonic (3{t). The result was one of the first maps on the plane of the 
dimensionless amplitudes (/3o,/3i), which differs notably from the Strutt diagramme 
(see Fig. |4|). 

By looking for Tr6(l) = ±2cos^, the computer identified also a sequence of 
curves in the stability areas where 



granting that oscillatory motions close after the n periods of H{t), forming the twin 
loops in both 2D subspaces {x,Px) and {y,py). In order to generate the loop of 4 
canonical variables, one must assure that the rotation r(l) simultaneously closes. The 
harmonic component /3i sin(27rt) does not contribute to r{l). The only condition is 
that the constant intensity /Sq should rotate the canonical variables by ±2A;7r after 



X± = e 



(56) 
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Figure 4: The types of Floquet operations 17(1) generated by the harmonic case of 
(|47p on the (/3o,/3i) plane (the map of Delgado). The clear areas are the stability 
zones, in which the evolution loops can occur for Trb(l) = ±2cos^, under the 
subsidiary condition r(l) — 1. 



some m repetitions {k,m = 1,2,...). This distinguishes the sequence of straight lines 
Po = on the stability map of Fig S) Their intersections with the loop curves 

yield the amplitude pairs (/3o,/3i) generating the loop phenomena for all 4 canonical 
variables after a finite number of mn repetitions (see Fig. [5]). 

While the existence of the harmonic loops is known [89|, their extremely regular, 
kaleidoscopic forms have some more implications. One of them is the exact vanishing 
of the "operator centers" X = {X,Y). Indeed: 

Proposition 2 Suppose, for a pair of amplitudes (/3o,/3i) the evolution loop of the 
oscillatory part b(t) closes for t = n, while the rotation r{T) yields a certain non- 
trivial angle 9 = (^) 2tt, where ^ is rational but not an integer (no matter, whether 
^ is smaller or greater than 1). Then, the loop obtained for t ~ mr ~ nm by 
simultaneous closing of both, oscillatory and rotational motions has the trivial center 
on both classical and quantum levels. 

Proof. In fact, the rotation r(r) breaks the oscillatory loop at t — t, marking a 
new end point qe rotated by 9 with respect to the initial q. The qe , in turn, becomes 
the initial point for the next fragment of the trajectory, which again does not close, 
but is just the ^-rotated version of the previous one, ending up at q2e, etc. As the 
result, the whole trajectory is the sum of similar fragments, generated between the 
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Figure 5: The classical trajectories in the xy plane illustrating the evolution loops. 
a) For P{t) — — J — 1.152 sin {2nt), closing after 15 periods. In blue the deformation 
(but not drifting) under the constant force F — ( — 1.5, 0). b) The evolution loop for 
/3(t) = f — 0.815 sin {2nt) closing after 24 field periods. In green, the deformation 
by the constant force F = (1,0). 



time moments 0, r, 2r, ... joining the subsequent points q, qe, q2e, • • ■, each one just the 
^-rotated version of the previous one, until finally, at t = mr, the trajectory closes 
with q^e = q- As the sum of the 0-rotated steps, the whole loop is invariant under the 
0- rotation and so is the loop center X defined by ([25l) . However, the only vector (with 
numerical or operator components) invariant under a non-trivial rotation around the 
coordinate center is X = O.B 

Consistently with Sec. [5] this assures the stability of the harmonic loops which 
resist drifting (see Fig. [5]) . It looks almost as an elementary analogue of the anomalous 
resistance observed in 2D [111] (though the attempts of overestimating the analogy 
might be misleading). 

7.2. The biharmonic case 

In turn, we shall consider the biharmonic fields (pS)) or (|T7|) with /?o = 0, defined 
by pairs of the dimensionless amplitudes /3i , /32 ■ Their agreeable property is that 
at each t = nT = n all rotations cancel, and simultaneously, the field amplitudes 
I3{n) = 0. Henceforth, the evolution operators U{nT) = U{n) are reduced to the 
purely oscillatory part; the corresponding one-period evolution matrix (|5ip reduces to 
a pair of &(1) cells. Moreover, the physical sense of the corresponding evolution in the 
time moments nT = n (n e Z) is not affected by the difference between the canonical 
and kinetic momenta. 

The stability areas determined by the computer scanning of Trfo(l) (c.f. |112) ) 
form now a new 2-parameter map (see Fig. [71). As before, they are densely populated 
by curves Tr6(l) — ±2cos(^^), whose points (/3i,/32) generate loops, though now 
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Figure 6: Biharmonic evolution loop for /3(f) — |- sin(27rf) + 9.966 sin(47rf) on 
x,y-p\a,ne. The loop closes after 6 field periods. 



(/3i,/32) don't need to fulfill any extra condition as the rotations r(l) automatically 
cancel for t = n E Z. A loop which closes up after 6 field periods is drawn on Fig. [6] 
Note its symmetry under the parity reflection, which grants the exact vanishing of X. 
However the parity is not even necessary. 

Proposition 3 Consider a periodic sequence of field pulses (3(t) = (3(t + 1), with 
j3{k) = 0, for k <E J., and moreover, f3(t)dt = (i.e., the rotations cancel in each 
full period T = 1 and its multiples) . Then the evolution loops which close for the first 
time after a certain number n > I of the (3 -periods has the vanishing center X = 0. 

Proof. Within our assumptions the eigenvalues X± of 6(1) are the n-th roots of 
unity, A" = A" = 1. If now A-|- ^ A_, then 6(1) is diagonalizable, with complex 
eigenvalues, and so is the 4x4 matrix u{l). If A+ = A_ the only possibilities are 
A-i- = A_ = 1 and A-|- = A_ = —1. Should then b{l) be non-diagonalizable, it could 
not produce the closed process, contrary to our assumption. Hence, the cell must be 
trivial and 6(1) diagonalizable. If so, A+ = A_ = 1 6(1) = 1 is excluded by the 
assumption that the loop closes only after a number n > 1 of steps. Hence, either 
6(1) = —1, but if not, then 6(1) is a real matrix with a pair of different complex 
eigenvalues. The corresponding u{l) has either 4 identical eigenvalues (—1) or two 
(identical) essentially complex pairs. The evolution matrix u{t) in the sequence of 
time intervals [0, 1) U [1, 2) U • • • U [n — 1, n] reproduces the u{t) from [0, 1), but pre- 
ceded by the increasing powers of w(l): 



u{l) 
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If now m(1)" = 1, then the whole matrix trajectory is invariant under the muhi- 
phcation by u(l) and so is the trajectory center X which must fulfiU u{l)X = IX; 
but since u{l) has no eigenvalue 1, the center must vanish identically, implying the 
loop stability (the loop resists the external forces). ■ 

The group theoretical sense of the 2D-loops can be noticed by reinterpreting the 
time variable 9 = 2'Kt as the rotation angle in the plane. The quantity which returns 
to its initial value after a full 27r-rotation is a tensor. An entity which changes its sign 
is a spinor. An entity which returns to its initial value only after a finite number of 
27r-rotations (as e.g. the loops of Figs. [5] and |6]) obeys a fractional representation of 
the rotation group in 2D |113j . 



7.3. The imperfections 

To estimate the corrections to the scheme it is convenient to return to the ordinary 
units. For slow classical/quantum motions (low kinetic momenta), the use of the 
Schrodinger's QM seems justified; but some dissipative correction must affect our 
semiclassical results. If the pulsating field is exactly periodic, the probability of 
quantum jumps might be read from the quasienergies of the Floquet Hamiltonian. 
Yet, for the low frequency of the applied fields and the evolution limited to a few pulse 
periods, more reliable seems the Abraham-Lorentz radiative acceleration arad = cr q , 
where a = is the "characteristic time" of the charged particle (c.f. Jackson 

|114j ). Taking as the model the electron with a = 2.12 x 10~^'*sec one has in the 
dimensionless variables: Qrad = y 9 ; so in the harmonic field frequency corresponding 
to the long w 1km RW the single operation period T = ^ 3.3 x 10"^ implying 
very little A-L corrections arad ~ 6.4 x 10^^^ Q (though the evaluation changes for 
increasing uj; see also [96]). 

A different source of errors are the retardation free vector potentials in the 
laboratory approximations (1361) . Their relativistic equivalents can be obtained by 
looking for the freely propagating fields of cylindrical symmetry; 



A„(t,r) = sinwt (I)(r) 



X 



(57) 



rO + 43> = (58) 
dr 



The d'Alembert eq. DA^ = then traduces into: 

(compare [115] ) which admits an analytic solution (S>{r) — 00 + 02 ^^ + 04 r* + - • • with 
the coefficients given by the recurrence 02„ = - (|j) „(„+!) ^2(k-i) . n = 1,2, • • •. 
They show that the "laboratory Hamiltonian" (j44M8|) is once again the first step of 
the EIH method [53] for the potential (|57II58I) . Its precision depends on the ratio 
between the effective size Ar of each loop and the wavelength associated with the 
pulse frequency uj, i.e. A = cT = c^. If A is much greater than the loop size, 
^ 1, then the formula (|57]) describes the harmonic magnetic fields differing little 
from the homogeneous field (|46|) , drawing with a good accuracy the regular loop forms 
predicted by the quadratic Hamiltonian 1^^. Note that an analogous estimation (of 
the ratio between the trajectory size and the field inhomogeneity parameter) works 
also in different circumstances, e.g. for the atomic or molecular systems irradiated by 
the coherent laser light. The typical atomic size {e.g. of widely studied Rb atoms) 
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is in the range w 10~^cm(~ 1^) but the A of visible Ught is much greater (between 
3800 A and 7800 A); which seems one of the reasons why the approximate description 
of the electric forces by a homogeneous, pulsating field E ^ Eq sin (out) gives so good 
results in the description of Rabi rotations |116) . The question whether the effects 
we described can have some analogues for the exact relativistic solution in strongly 
inhomogeneous fields (e.g. |117) ) is open. 

Until now, our exploration was limited to the stability areas with their typical 
"fauna" and "flora" of drifting and stable loops. However, the problems of quantum 
control unavoidably leads us across the borders. 



8. The threshold: (5-kicks and inverted free evolution. 



Since the rotations commute with the whole Hamiltonian (|48p. we shall again fix 
attention on the b{t) cells of the evolution process (|53l) . 

One of the major control challenges are the idealized (5-kicks of the elastic 
potentials and the events of the "distorted free evolution", both represented by 
extremely simple operators and (7,p-transformations represented by 6(1): 



1 



1 


T 





1 



(r^T), (59) 



(where q, p mean x^p^ or y,Py) though almost impossible to achieve straightforwardly. 

The left operation in (j59p can be interpreted as a formal result of an infinitely 
strong pulse of an oscillator potential applied within an infinitesimal time. The right 
one, if r > 0, r ^ 1, represents the distorted (slowed or accelerated) free evolution; 
but if T < 0, it yields highly counter-intuitive effect: it inverts the free evolution of a 
wave packet, including the interference between the different packet components. No 
less provocative are the evolution operators and matrices: 

1 



Pe- 



1 



Pe- 



1 


T 




1 



(60) 



i.e. ([5^ superposed with the parity reflection P : q ^ —q,P — >■ —p- AH matrices ([5^ . 
(|60p belong to the resonance borders (case II of our classification). In principle, they 
can be generated by combinatorial field patterns (TT] |31] 1118) . However could they be 
caused by a "soft persuasion" instead of brutal force? 

With some surprise we found that (|59p and/or (j60p occur on the separatrix 
branches of Fig. [7] with either Tr 6(1) = +2 (we call them of type or Tr 6(1) = —2 
(type — ). The result contained some mystery. Though it is quite obvious that for 
|Tr6(l)| = 2, the matrix 6(1) degenerates, this still does not imply that it must appear 
in the standard Jordan's form ([59|l or (|60|) precisely in the q,p basis. 



As turns out, the phenomenon occurs always when /3(i)^ in eq. (|52p is symmetric 
with respect to the center of the operation interval. To show this it is convenient to fix 
the time variable t with < = at the symmetry center; then to consider the evolution 
matrix b{t) ~ b{t, — t) generated in the expanding interval [— (see [96]). 

Proposition 4 Suppose A(t) is symmetric around the point t = 0, i.e. A(— t) = A(t). 
Then if at any t G R, Tr b{t) — 2, the matrix b{t) = b{t, —t) must take one of the forms 
I159\) . but if Tr b{t) — —2, then b{t) must take one of the forms i60\) . 

Proof. For the symmetric A{t) given by (|5^ the evolution matrix 6 ~ b{t, —t) fulfills: 
^ = A(i)6 + 6A(i) = (621 - /3'6i2)l + (Tr 6) A(i) (61) 




Figure 7: The map of the stabihty, resonance areas and the separatrix borders of types (+) and (— ) on the intensity map /?i,/32, scanned by 
integrating the matrix equation (|52p for the biharmonic 6(1) in the periodicity interval [0,1]. The separatrix borderlines host the exceptional 
matrices of types (|59p . (|60p . The branches yielding the oscillator kicks are brown, and the ones representing the distorted free evolution are green. 
The approximate values of the distorted evolution time r and the simulated oscillator kicks a are given by the sequences of numbers along the 
corresponding borderlines. Notice the negative r on two sections of the internal green branch. The points A, B, marked blue and green represent 
the field amplitudes generating the space trajectories of Fig. 6 and Fig. 8 respectively. 
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Hence one easily shows: 

= ^ 612621 - ^(Tr 6)2 = C = const. (62) 

The constant in (j62|) can be determined by shrinking the interval [—t, t] to zero. For 
t = 0, the evolution matrix is trivial, 6 = 6(0,0) = 1,612 = 621 = 0, i.e., the constant 
in (|62]) is C = -i(Tr6)2 ^ -1, implying 

612621 = ^(Tr 6)2-1. (63) 

Hence, 612621 must vanish whenever Tr6 = ±2. So, either 612 or 621 must vanish, 
leading to one of Jordan's forms (ICT)) or (|60p on borders of type (+) or (— ).■ 

Corollary. While under our symmetry assumptions the values 612 — and 621 — 
appear on the stability thresholds, none of them can occur inside of the squeezing ar- 
eas. Indeed, (|63p imply that for \Tr 6| > 2 none of 612 and 621 can vanish. 

Note that the Proposition 4 applies to the biharmonic (3{t), odd with respect 
to the centers of the periodicity intervals [n,n + 1], and vanishing at the extremes 
n = 0, ±1, ±2, • • [jj. The corresponding stability borders (Fig. [7|) form a sequence 
of closed, intersecting pairs of branches surrounding the map center. The computer 
scanning shows that the inner pair, intersecting at 4 points, is of type (— ), i.e., hosts 
the evolution matrices (j60p . The next pair, with 2 intersections, is of type (+), hosting 
(|59p . In general, the 1**, 3'"'*, and all odd pairs are of type (— ), whereas all even pairs 
are of type (+) (see Fig. [7]). On each pair, one branch (brown) offers the soft imitations 
of the (5-kicks (with the amplitudes a 7^ 0, but t — 0), while the other one (green), 
yields the incidents of the distorted free evolution {a — 0, t ^ 1). Note, that the 
operations on the (+) branches give the clean effects after the single pulse period (or 
its multiples), while the operations of the (— ) branches yield the parity free effects 
2a or 2t after two pulse periods (or 2na, 2nr after 2n repetitions). To illustrate the 
details, our Fig. [7]reports the values of the "distorted time" t along the green branches 
and the amplitudes of the effective oscillator kicks along the first brown branch. 

Our computer scanning shows that the double period operations of the type (— ) 
have the exactly vanishing Floquet centers, and so, produce the stable (resistant) 
effects, while the operations of type (+) define non-vanishing fuzzy centers and can 
be affected by constant external forces. 

As an example, we have choosen a case of the free evolution inversion for a point 
on the first (negative) separatrix branch of Fig. [T] The effects of the "retrospective 
operations" after 4 double periods are illustrated by the center of a wave packet 
recovering in the time instants t = 2n (n = 1, 2, • • •) its past shapes (shrinking instead 
of spreading in some number of steps); its subsequent positions showing a sequence of 
shifts in the direction opposite to the initial velocity (Fig. |8]). 

In turn, the possibility of simulating the S{t) pulses of the oscillator forces by 
the soft biharmonic fields (brown branches of Fig. [7]) is no less essential: it saves the 
realistic sense of all operations programmed with the help of the sharp oscillator pulses 
[72l|80l|73]. 

II For the harmonic 13 (t), the assumptions of Proposition U] are satisfied in the periodicity intervals 
[n — i,n -I- i] thanks to the symmetry of (3{t). However, due to the existence of the constant field 
component, the vector potential Ave does not vanish at the operation extremes, affecting the physical 
sense of II59I I and II6OI I (c.f. the remarks in Dodonov et al. |31|). 



dt 



612621 



^(Tr6)2 
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Figure 8: The double period operation for tlie biliarmonic field amplitudes 
(/Si, (82) = (2.40,2.68) on the first separatrix branch. Under the repeated 
applications of the 2 period pulse pattern, the Gaussian packet shrinks instead 
of expanding and its center travels in the direction opposite to its initial velocity, 
simulating the incidents of the inverted free evolution with At = nr ^ — 0.3688n 
at n = 2,4, 6 - • •. 



9. Dark areas: the squeezing 

Below we shall not comment any particular squeezed states, but rather the squeezing 
operator which can deform any wave packet, no matter, is it Gaussian or not, and we 
shall be looking for the origin of this phenomenon in the structure of the full period 
evolution (Floquet) operator ?7(1) and its matrix u{l). 

The squeezing mechanisms were carefully examined in an ample sequence of pa- 
pers, the most complete ones by Dodonov etal[2Il[SIllSllSllMllSllS3II3ISl ITOTl 
I1051I1121I1181I119) : so our approach cannot dissent too much, though it can contribute 
with some more observations: 

The use of the quadratic Hamiltonians with the time dependent B = B{t) is justi- 
fied as the first step of the known EIH method [93] for not too ample experiment areas. 

In case of periodic /3(i) = /3(i -t- 1) in the cylindrical geometry the information 
about the squeezing can be read from the 2x2 oscillatory matrix 6(1), even if the 
time dependent field pulses do not distinguish any invariant "squeezing center" . This 
holds also if the magnetic field vanishes at the beginning and at the end of the opera- 
tion. If |Tr6(l)| > 2 (dark areas), then A± are real, A+A_ = 1, describing a squeezing 
on the canonical plane. For Tr6(l) > 2 (type +), the transformation is a squeezing 
sensu stricto with two positive coefficients X± > 0, but if Tr&(l) < —2 (type -), then 
6(1) yields the superposition of squeezing and the parity transformation (see also our 
recent report [96'). The squeezed or amplified canonical observables are defined just 
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by the (row) eigenvectors of 6(1) (c.f. 

Some studies dedicate special attention to the scale squeezing or expansion in 
which q is squeezed at the cost of p or vice versa [Ml I118| . Here, the cases of 612 — 
or 621 =0 are quite relevant not only on the separatrix branches, but also inside of 
the squeezing areas. The value &12 = in 6(1) implies the coordinate q squeezed or 
expanded at the cost of some other canonical observable, while 621 = 0, the same for 
the canonical momentum. The coincidence 612 — 621 — 0, would represent the position 
q (meaning x and/or y) squeezed or amplified exactly at the cost of the corresponding 
momentum (p^ or py). However, due to our Proposition 4, none of these cases can 
happen if is symmetric around the operation center, when the desired zeroes of 612 
and/or 621 are all located on the separatrix, without penetrating ever inside the dark 
areas. Hence, the scale transformations can be achieved only in the operation intervals 
with an asymmetric /3(t)^. In this last case the role of the 'fuzzy points' again calls 
attention. If the Floquet center X of the classical/quantum trajectory vanishes, then 
the squeezing operations will be immune to the constant external force F. However, 
should X be fuzzy, then if the field pulses are repeated, the g-trajectory, apart of the 
sequence of the squeezing operations, will show a drift in a direction orthogonal to 
F. So, would the "fuzzy center" be indeed a generalized element of noncommutative 
geometry applicable for variable fields? 

10. Open Problems. 

Since the effects which we discussed are computed in the dimensionless field variables, 
it is not excluded that their analogues can reappear in different orders of magnitudes 
as e.g. interiors of the oscillating cristal lattices, the nodal points of the crossed 
electromagnetic waves |1201 [80] . which do not involve the macroscopic traps |121) . 
Supposing that some analogues 'in little' can exist, they might contribute to different 
topics recently studied. 

The behavior of particle gas in 2D in presence of the harmonic or biharmonic 
fields, shows a sequence of surprising reactions to the external electric force, in form 
of a "giant" |122| or vanishing, or even negative resistance |123) . The effects are 
apparently due to the presence of material shells forming a sandwich on both sides 
of the motion plane; but the theoretical discussions are not yet concluded (see e.g. 
Zudov et al. [111) ). One might henceforth ask, whether the surprising phenomena 
have not some elementary counterparts for time dependent fields? Could e.g., the 
behavior of the 2D evolution loops with the vanishing 'fuzzy center', be interpreted 
as an analogue of an 'anomalous resistance' ? Could the distorted free evolution mean 
a modified effective time? Could the squeezing effects matter? Last not least, could 
the analogous phenomena affect the polarized vacuum, if such an entity exists |124j ? 
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